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Objectives 



Entry Systems and Technology Division 


Primary 

To predict surface distributions of pressure and heat flux using “standard” 
simulation model(s) for: 

(a) Sharp 25° /55° double cone model 

(b) Hollow cylinder-flare (30° ) model 

tested at laminar flow conditions in LENS-XX at CUBRC 


Secondary 

To understand/quantify sensitivity of predictions to uncertainties in freestream 
conditions 

To understand/quantify the influence of physical models (thermochemical vs. 
chemical nonequilibrium) on flow predictions 


Focus of this presentation is solely on the double cone model 
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Double Cone Model 



Entry Systems and Technology Division 


o 



Cone-flare model has a sharp tip 
Same model has been tested in LENS-1 (reflected shock tunnel) 
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Modeling and Simulation 




Simulation Tool 


V4.03.1 of Dplr (axisymmetric) 


Physics 


•5-species air (N 2j 0 2j NO, N, O) for all cases 
• 2 Physical models 

- Thermochemical nonequilibrium - TCNEQ 

- Chemical nonequilibrium - TEQ 
•Thermodynamic properties of species: Lewis curve fits 

- For consistency with Cheetah 

- Typically would have used SHO/RR 

•Momentum and energy transport properties: Gupta-Yos curve fits 
•Mass transport: Self-consistent effective binary diffusion 
•Chemistry: Park 90 mechanism and kinetics 
•Wall model: Fully catalytic to atom recombination and T w = 300 K 


•1 st -order accuracy in time and 2 nd -order accuracy in space 


Numerics 


Grid tailoring capability of Dplr exercised in all computations 
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Strategy 



Entry Systems and Technology Division 


• Establish grid requirements via the following studies: 

- Wall-normal grid resolution 

• Cell Reynolds number used the criterion 

• Final cell Re used is 0.1 

- Wall temperature sensitivity 

• Heating results insensitive to choice of wall temperature for high total 
enthalpy cases 

• At lowest total enthalpy, less than 1% change for wall temperatures beween 
400 and 300 K 

- Grid convergence established by computing on grids that are 4x coarser and 
finer than nominal grid 

• Results reported on nominal grid only 

- 2 nd -order spatial accuracy found to be adequate for all axisymmetric 
computations 

• Some sensitivity seen in 3D computations, but 3D grids employed were 
coarser than nominal axisymmetric grid 
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Freestream Velocity Uncertainty Line of Inquiry, I 


Entry Systems and Technology Division 


• Surface heating is proportional to the “enthalpy potential” 

- q JLX H 0 - h w 

• Total enthalpy is a conserved quantity across a normal shock 

- Inviscid flow is adiabatic 



Total enthalpy is a sum of flow kinetic energy and thermochemical enthalpy 

ex p( JT,) 


H 0 (T,T v ) = ±V 2 + X 


s=molecules 


RT 

3.5 + 

f \ 
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s=atoms 


If freestream is undissociated air, then no atomic contribution 

ex p( w , /r ) 


H 0 (T,T V ) = \v 2 + 


c s RT 3.5 + v ’* 


5— N 2 ,0 2 


T 1 ex P( v,,_ IT,. ) 


Freestream kinetic energy is the dominant component - Velocity is the sensitivity variable 
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Freestream Velocity Uncertainty Line of Inquiry, II 



Entry Systems and Technology Division 


• Pitot pressure can be measured reasonably well - P 0 = p TO \fj 

- Cannot break this product without an independent measurement of velocity 

• Assume uncertainty in given freestream velocity is ±x%, but P 0 is fixed 


j^mew 


old 


f = l±x 


Therefore, 


pnew pold 

M) ~ r Q 


^ ynew ^ “ old ^ -yold ^ ' 


1 


old 


f 


Still have to specify T* 

- One way is to assume static pressure is also fixed (at nom.) 


„ new old 

P =P 


new ^ . rj-mew old yyyold 

s = 1 


rj-mew ryold 


Freestream enthalpy will change! => Freestream velocity and total enthalpy are synonymous 


H 0 (r ew ) = f -{V oU f + 3.5RT oU =f 2 H„(T° u ) 
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Learning Cases - Cases from LENS-1 Matrix 

(“Open” Validation Cases in AIAA 2013-2836) 


Entry Systems and Technology Division 


Run 91 Run 8 


Ao/g-nr 3 

1 .8342 

1 .5498 

1.291 

1.34 

VJkm.s' 1 

2.731 

4.148 

3.067 

4.063 

tjk 

190.1 

729 

166 

303 

tjk 

1001 

773.1 

2711 

3085 

Gas comp. 

(mass fractions) 

0 2 = 0.9986 
O =0.0014 

0 2 = 0.9389 
O =0.0611 

N 2 = 0.9999 
N =0.0001 

N 2 = 0.9973 
N = 0.0027 

Ho/MJ.kg 1 

3.99 

10.26 

5.28 

9.17 

Molecular (tr/rot) 

4.3% 

6.1% 

3.3% 

3.4% 

Molecular(vib/el) 

1 .7% 

0.3% 

7.6% 

5.5% 

Atomic 

0.5% 

9.7% 

| 0.1% 

1 .0% 

Kin. energy 

93.5% 

83.9% 

89.0% 

90.1% 



Freestream conditions from numerical simulations of nozzle flow 
Double cone configuration has been extensively studied by Nompelis et al. (AIAA 2010-1283) 
Unable to obtain a stable solution for Run 80 with TCNEQ model, but TEQ model is stable 
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Learning Case - Run 90 (02/0 Mixture) 



Entry Systems and Technology Division 


Pressure 


Heat Flux 



Predicted environments relatively insensitive to choice of TCNEQ or TEQ 
Velocity has a bigger impact - lower velocity larger bubble size 
Velocity influence on the 55° cone is also large 
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Learning Case - Run 90 (02/0 Mixture) 



Entry Systems and Technology Division 


Pressure 


Heat Flux 



If velocity is 95% of nominal, then better replication of pressure data 
Error bars on experimental data not provided 
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Learning Case - Run 91 (02/0 Mixture) 




Entry Systems and Technology Division 


Pressure 


Heat Flux 


o 


o 

Experiment (Run 91) 


V -1.0V nom (TEQ) 


- V =1.0V nom (TCNEQ) 


V = 0.9 V 




o 

LO 



0 20 40 60 80 100 120 140 160 

Axial coordinate/mm 


o 



If velocity is 90% of nominal, then better replication of pressure data 
The double peak characteristic post-impingement needs schlieren image for confirmation 
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Learning Case - Run 91 (02/0 Mixture) 



NASA 


Entry Systems and Technology Division 


Nominal Velocity 


90% Nominal Velocity 




Velocity magnitude determines whether or not there is a Mach stem 
A schlieren image for this case would be most helpful 
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From “Open” to “Blind” Validation Cases 



Entry Systems and Technology Division 


• “Open” validation cases are from LENS-1 

- LENS-1 is a reflected shock tunnel 

• “Blind” validation cases are from LENS-XX 

- LENS-XX is an expansion tunnel 


Still retain the freestream velocity uncertainty line of inquiry 
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“Blind” Validation Cases - LENS-XX Matrix 




Entry Systems and Technology Division 



Case 1 

Case 2 

Case 6 

Case 5 

Case 3 

Case 4 

Ao/g-nr 3 

0.499 

0.984 

2.045 

1.057 

0.51 

0.964 

1/Jkm.s- 1 

3.246 

4.303 

5.466 

5.996 

6.028 

6.497 

tjk 

175 

389 

573 

523 

521 

652 

P 0 / kPa 

5.1 

17.5 

59 

36.8 

18 

39.5 

Ho/MJ.kg 1 

5.44 

9.65 

15.23 

18.51 

18.7 

21.77 

Thermal 

3.2% 

4.1% 

3.7% 

2.8% 

2.8% 

3.0% 

Kin. energy 

96.8% 

95.9% 

96.3% 

97.2% 

97.2% 

97.0% 


12.2 

10.9 

11.46 

13.14 

13.23 

12.82 

Re uoo xio- 

6 /m _1 

0.14 

0.19 

0.39 

0.23 

0.11 

0.2 

Aoo/mm 

0.129 

0.085 

0.043 

0.084 

0.178 

0.095 

5 

8 

X 

O 

CO 

1.40 

0.92 

0.47 

0.91 

1.93 

1.03 


Freestream conditions from CHEETAh code (OD/Equilibrium) 
Code calibrated to shock velocity and pressure measurements 
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Blind Study - Case 1 (Low Enthalpy) 

Influence of Freestream Velocity 



Entry Systems and Technology Division 


Pressure Heat Flux 




The test configuration is shown in light grey lines 
from "keep M fixed" or "keep Re fixed" idea has little influence on predictions 
±10% variation in velocity at low enthalpy (5 MJ/kg) has more influence 
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Blind Study - Case 4 (High Enthalpy) 

Influence of Freestream Velocity 



Entry Systems and Technology Division 


Pressure Heat Flux 




from "keep M fixed" or "keep Re fixed" idea has little influence on predictions 
±10% variation in velocity at high enthalpy (21.8 MJ/kg) has more influence 
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Blind Study - Cases 1 & 4 

Influence of Angle of Attack 



Entry Systems and Technology Division 



Low Enthalpy 
(Case 1) 


a = 0° 



0 15 30 45 60 75 90 105 120 135 150 



a = 2' 


High Enthalpy 
(Case 4) 


a = 0° 



q/W.cm ' 2 


200 400 600 800 1 000 1 200 1 400 1 600 1 800 2000 


Angle of attack influence is greatest for largest separation bubble size 
Bubble size at 5 MJ/kg » Bubble size at 21.8 MJ/kg 
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Blind Study - Case 1 (Low Enthalpy) 

Influence of Gas Model 



Entry Systems and Technology Division 


Pressure Heat Flux 




SHO/RR or LeRC fits for thermodynamic properties makes no difference 
Thermal relaxation at low enthalpy (5 MJ/kg) has an enormous influence 
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Blind Study - Case 4 (High Enthalpy) 

Influence of Gas Model 



Entry Systems and Technology Division 


Pressure 


Heat Flux 



o 

o 



SHO/RR or LeRC fits for thermodynamic properties makes no difference 
Thermal relaxation at high enthalpy (21.8 MJ/kg) has smaller influence than in Case 1 
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Concluding Remarks 




• Accomplishments 

- All cases computed for both configurations 

• The hollow cylinder-flare configuration has not received the same level of 
attention as the double cone configuration 

- Each case had something unique to offer, especially Case 2 


• Things still left to do 

- Study influence of transport properties 

- Complete grid convergence studies for the hollow cylinder-flare configuration 

• Lesson Learned and/or Open issues (in the view of the author) 

- Since freestream velocity has the most influence on environments, would be 
preferable to have independent measurements of freestream velocity 

• Would also be preferable to work with the actual nozzle contour and perform 
3D flow computations -> predicted freestream chemical state is a bonus 

- Need some cases from facility characterization as well to complete V&V story 


• Schlieren images for these cases would be most helpful 


• Stag, measurements of heat flux and pitot pressure for either the 3.5-in 
cylinder or 1.25-in hemisphere 



